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Abstract 

We use high precision lattice simulations to calculate the plaquette expectation value 
in three-dimensional SU(A'^c) gauge theory for — 2,3,4,5,8. Using these results, 
we study the iVc-dependence of the first non-perturbative coefficient in the weak-coupling 
expansion of hot QCD. We demonstrate that, in the limit of large Nc, the functional form 
of the plaquette expectation value with ultraviolet divergences subtracted is 15.9(2) — 
44(2)/iV2. 



1. Introduction 

The determination of QCD pressure up to order is a long-standing problem in finite- 
temperature field theory [B El E]- This is the first order where a coefficient of the weak- 
couphng expansion, due to infrared divergences, gets contributions from an infinite number 
of loop-diagrams and thus is non-perturbative. 

However, at high enough temperatures (T > 2Tc) the properties of finite-temperature 
QCD can be described by dimensionally reduced effective field theory methods [H [5]. By 
integrating out temporal degrees of freedom a three-dimensional pure gauge theory, called 
magnetostatic QCD (MQCD), is constructed. This allows us to isolate all the divergences 
to MQCD and study it using lattice calculations. The integration out is most conveniently 
performed perturbatively in MS scheme |6]. 

We can relate any lattice regularized quantities within MQCD to the continuum scheme 
(MS), because MQCD is super-renormalizable. There are ultraviolet divergences up to 4- loop 
level only j7j . Terms required in the conversion have been determined up to 3-loop level [U |9] . 
Infrared divergences cause an additional complication in the 4-loop level. The computation 
requires an introduction of an IR cutoff, which then cancels once lattice and MS results are 
subtracted. This computation has been carried out recently for A'^c = 3 in [10] using stochastic 
perturbation theory. 

In [11] the plaquette expectation value, which determines the non-perturbative contribu- 
tion, was measured for Nc = 3. The purpose of this paper is to extend the results to study the 
A'c-dependence of this observable. We carry out lattice measurements of the plaquette with 
Nc = 2, 3, 4, 5 and 8 to obtain the A'c-dependence. We also get an independent approximation 
for the A'^c = 3 result. This acts as a consistency check for the whole pressure calculation. 
Namely, we expect to see smooth A'c-dependence in the observable. 

Additionally, there are various other physical motivations to study the A'c-dependence and 
especially the large- iVc limit of SU(A'c) gauge theories [12|. The limit ^ oo simplifies the 
theory significantly, but nevertheless the phenomenology is in many ways similar to SU(3). 
These reasons have motivated numerous large- A'c limit studies on the lattice [131 114j . 

The paper is organized as follows. In Sec. 2, we give the theoretical background of our 
study and specify the observable we consider. In Sec. 3 we present the numerical results of 
lattice Monte Carlo simulations. Conclusions are given in Sec. 4. 

2. Theoretical setup 

The ultimate interest of our study is Euclidean pure SU(A''c) Yang-Mills theory, defined in 
continuum dimensional regularization by 




(2.1) 
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where d = 3 - 2e, is the gauge couphng, k,l = 1, . . . ,d, Fki = i[Dk,Di], Dk = dk - iAk, 
Af^ = A'^T"', and T" are Hermitean generators of SU(A'^c) normahzed such that Tr[r"T^] = 
jl. The vacuum energy density in MS (suppressing Faddev-Popov and gauge fixing terms) 
is defined by 
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where V denotes the d-dimensional volume. The use of the MS dimensional regularization 
scheme removes any 1/e poles from the expression. In fact, using dimensional regularization 
the perturbative result vanishes, because there are no mass scales in the propagators and 
therefore the UV and IR divergences cancel each other. However, for dimensional reasons, 
the non-perturbative form of the free energy is 
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where p, is the MS renormalization scheme scale parameter. The coefficient of the logarithm 
has been calculated by introducing a mass scale uiq for gluon and ghost propagators and 
sending after the computation [3l [15] : 
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where dA = - 1. The non-perturbative constant part Bq, which is a function of the 
number of colors, is what one would ultimately like to determine. 

Using standard Wilson discretization, we can write the corresponding action on the lattice 

as 
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where Pki is the plaquette, a is the lattice spacing and /3 = 2Nc/{agl). Hence the continuum 
limit is taken by /3 — > oo. Analogously to MS, the free energy density is defined on the lattice 
as 
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Dimensionally, the vacuum energy density consists of terms of the form g'^^a^ ^. Thus, 
approaching the continuum limit, we can relate fa and as follows: 

A/ = fa-fus (2-7) 

= fin ^ + C[) + C2% + C-A + C,gl (in 4 + C^) + 0{gla). (2.8) 

Taking derivatives of Eq. ()2.7p with respect to g'^ and using 3d rotational and translational 
symmetries on the lattice, we obtain the relation [11] 
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The relations between c,- and C,- are 
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ci = Ci/3 C2 = -^C2 C3 = -^C's 

C4 = -mid C4 = (^4 - - - 2 ln(2iVc). (2.10) 
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The first follows from a straightforward 1-loop computation: 

ci = f . (2.11) 
The 2-loop constant has been computed in three dimensions in [16] and can be written as 

C2 = — ^ + vrS + 4K1 + -K5 (2.12 

3 (47r)2 V3iV2 ^4 2 3 7 

= dAiV^ 1^0.03327444(8)-^^), (2.13) 

where the coefficients S, ki and K5 can be found in [71 [T7]. The 3-loop term has been 
computed in three dimensions recently in Ref. [9]: 

C3 = cLaN^ (^0.0147397(3) - 0.04289464(7)^ + 0.04978944(1)-^) . (2.14) 
Because there is no Jl dependence in fa, the value of C4 is determined by fy^, 

C4 = 0.000502301323dAiVc • (2.15) 

The four-loop free energy itself is an IR divergent quantity at in both MS and lattice 
schemes. But the finite difference between them, C4, can be defined by introducing the same 
IR cutoff, e.g. a gluon mass, to both schemes. The cutoff dependence then cancels out when 
the two schemes are compared. At present C4 is known only for A'c = 3, for which it has been 
calculated using stochastic perturbation theory [TOj. 

For later use we define the quantity 
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which is a normalized plaquette expectation value minus all the ultraviolet divergences. 
Hence, 

^G(iVc) - - ^vr^) 4 = PG(cx),iVc). (2.17) 

Our goal here is to determine Pg(oo, Nc). After the A'^c-dependence of C4 has been determined 
by, e.g., stochastic perturbation theory, one has reached the final goal, the determination of 
Bg{Nc). 
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Figure 1: The significance loss due to the subtraction of ultraviolet divergences in the pla- 
quette expectation value with different Nc- Here "plaq"= (1 — ■j^Tr[P]) and the symbols Cj 
in curly brackets represent which subtractions of Eq. (j2.9p have been taken into account. 



3. Lattice computations 

The simulations were performed using Kennedy-Pendleton quasi heat bath (HB) [18] and 
overrelaxation (OR) algorithms. For the overrelaxation we used an algorithm which updates 
the whole matrix using singular value decomposition and performs very well for large A^c jl9] - 
Lattices of size N^, iV = 24,... ,400 were used. 

For each HB update we performed one OR. The number of updated subgroups in HB for 
Nc = 3, 4, 5 and 8 were 3, 4, 8 and 24, respectively. These subgroups were chosen randomly 
for each update. After each of these cycles we measured the value of the plaquette. The 
integrated autocorrelation times were around 0.75. For SU(2) we used dedicated OR and HB 
algorithms, with a ratio of one OR step for each HB update. The autocorrelation time was 
around 0.6. The data sets used for SU(3) are the same as in [TT] . 

The contribution of Bq to the plaquette expectation value in Eq. (j2.9p is about five orders 
of magnitude smaller than the leading order contribution. Thus we experience massive signif- 
icance loss in the subtraction and the accuracy requirement makes the numerical computation 
demanding (Fig. [T]). 

The only physical scale in this problem is the correlation length of the lightest glueball, 
which according to [20] is ~ l/NcQ^- The requirement, that this scale be in the reach of the 
lattice gives us the condition 

a<^T— <iVa, (3.1) 

glNc 

which translates into 

2iV2 < /3 < 2NIn. (3.2) 

Systematic errors due to the finite- volume effects turn out to be well under control. Because 
the theory is confining, we expect finite-volume effects to be exponentially suppressed when 
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Figure 2: Pg{/3,Nc) as a function of the physical lattice size (3/{NN^). Points denoted 
by open symbols are relatively low-statistics small volume simulations, included in order 
to illustrate the exponentially suppressed finite volume effects. These are omitted in the 
extrapolation. Finite-volume effects become visible when f3/{NN^) ^ 0.2. The points on the 
vertical axis indicate the infinite-volume estimate, obtained by fitting a constant to data in 
the range P/{NN^) < 0.1. 



the condition ()3.2p is fulfilled. As seen in Fig. [2l the finite- volume effects are no longer visible 
within our resolution when (3 < 0.2N^N. 

In Fig. [3] the effects arising from finite lattice spacing can be seen. We experience a 
qualitative change in the behavior of the plaquette expectation value at /3 ~ A'^^. The 
plaquette expectation value as a function of volume and lattice spacing a is consistent with 
the assumption of correlation lengths being ~ / (3. 

After numerous test runs we use in our simulations the requirement 

N^<13< A(Ae/3)2, (3.3) 

which is also the case in 

The continuum extrapolation is obtained by fitting a polynomial Pci^c) = (^1+^2//? +'^3/?^ 
to the infinite-volume extrapolated data in Fig. H] for each Ac separately. This functional 
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Figure 3: The solid line indicates the continuum extrapolation obtained by fitting a second 
order polynomial to the infinite-volume extrapolated data. Points denoted by lighter color 
are omitted. The bulk phase transition point is around -/V^//3 ~ 0.9. 

form describes data quite well. The x^/doi values for = 2,3,5 are excellent but slightly 
discouraging for A'c = 4, 8. The fitted values are show in Tabled! Using only statistical errors 
of the fitting parameters would underestimate the uncertainties of the continuum values, 
because the fit is dominated by points far from the continuum limit. Inclusion of higher 
order terms to the fitting function changes the continuum extrapolations by about one sigma. 
Therefore we expect that the 1-sigma error of the continuum extrapolated value is comparable 
to 2-sigma error of the fitting parameter di. 

At the leading order in N^, our measurements agree with the prediction of planar diagram 
theory with Pg{Nc,oo), approaching a constant (Fig. [5|). To study the next order contribu- 
tions we fit polynomials bi + b2/Nc, bi + b2/Nc + b^/N^ and bi + b^/N^ to the continuum 
extrapolated data in Fig. [6l We find that two last forms fit the data quite well. The 62 coef- 
ficient is zero (within our resolution) as could be expected from the form of the perturbative 
coefficient^, which are also functions of A'^. The data is not accurate enough to determine 
higher order terms. 

As our final results we quote 



""^Note, however, that terms ~ 1/A'c appear to be possible in certain other pure gauge theory observables [14] . 




= Pg(cx), N,) = 15.9(2) - U{2)/N^ (3.4) 



Inserting Ac = 3 we get 
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Figure 4: Continuum extrapolations of infinite- volume extrapolated data for each A^c- 





fit 


xVdof 


Pg{oo,N,) 


2 


5.09(15) 


+ 


16(3)/?-^ 


+ 




3(ll)r' 


5.1/6 


5.1(3) 


3 


10.7(2) 
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46(7)/?-! 


+ 


4.85(6) 


X 10^/3-2 


5.8/6 


10.7(4) 


4 


13.38(13) 


+ 


1.05(9) X lO^/?-! 
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2.58(14) 


X lO^p-"^ 


12.3/5 


13.4(3) 
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14.8(2) 
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1.8(2) X lO^/?-! 
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7.9(5) 


X lO^/?-^ 


7.7/4 


14.8(4) 
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14.7(2) 
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7.7(5) X 102/3-1 
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5.3(3) 


X 104/?-2 


17.7/4 


14.7(4) 



Table 1: The fitted values and x^/^^Qf of continuum extrapolations for each Nc- The value 
in the brackets indicates the uncertainty of the last digit. The last column indicates the 
continuum limit with systematic errors included. 
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Figure 5: Comparing the leading order behavior of Pg(oo, Nc) in Nc- As predicted by planar 
theory, Pq approaches a constant in the large- A'c limit. 



15 



10 

8 



5 







■■. ' 1 ' 1 









O Pq(oc,N^) 






b,+b,/N 

1 2 c 






— - b,+b,/N +b,/N ^ 

1 2 c 3 c 

— b,+b7N " 

1 3 c 













0.05 0.1 0.15 0.2 0.25 
1/N^ 

c 



Figure 6: Comparing different fits for higher order terms in Nc- The term is zero within 
our resolution implying that Pq is a function of N~'^. 
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which is consistent with the direct determination 10.7 ± 0.4 



4. Conclusions 

The purpose of this paper has been to measure the A^c-dependence of the expectation value 
of the plaquette in three-dimensional pure gauge theory. We have also outlined how the con- 
tinuum MS scheme free energy can be extracted from it. High precision lattice measurements 
of plaquette were performed with Nc = 2,3,4,5 and 8 and the large- A'^c limit was taken by 
extrapolation. We found that the non- perturbative input is Pq = 15.9(2) - U{2)/N^. The 
data does not seem to allow for terms ~ 1/-/Vc, and higher order terms, 0{1/N^) or 0{1/N^), 
are small enough such that the physical case A'^c = 3 is very well described by this form. 
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bi + b2N-^ 
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Table 2: Different fitting functions for Pg{oo, Nc). The term N^^ provides a very bad 
description of the data (1st case) or has a coefficient consistent with zero within our resolution 
(2nd case); see also Fig.[6j The confidence values of fits are plausible for the last two functions. 
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A. Tables 

In this Appendix we collect the numerical results for the plaquette expectation value mea- 
surements, which have been used in the continuum extrapolations. The column A'^ind gives 
the number of independent measurements within a data set. The data sets used for SU(3) 
are the same as in |llj . 
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SU(2) 




volume 


Nind 




6 


48^ 


40719 


0.1752161(16) 


7 


48^ 


42164 


0.1488698(13) 


9 


48-'^ 


43187 


0.1145493(10) 


9 


320^ 


5104 


0.11454906(17) 


11 


48-'^ 


42993 


0.0931322(8) 


11 


320^ 


8463 


0.09313207(11) 


13 


48-'^ 


44195 


0.0784776(7) 


13 


320^ 


8024 


0.07847755(9) 


16 


643 


157777 


0.06350205(19) 


16 


320^ 


14881 


0.06350198(5) 


20 


64^ 


271054 


0.05062861(11) 


20 


320^ 


12613 


0.05062829(5) 


24 


48^ 


904993 


0.04209730(8) 


24 


64^ 


317058 


0.04209720(9) 


24 


320^ 


14961 


0.04209733(4) 


32 


64^ 


868436 


0.03148821(4) 


32 


320^ 


15064 


0.03148828(2) 


32 


400^ 


6193 


0.03148828(3) 
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(0 

p 


volume 


A T 




12 


24^ 


13459 


0.2417125(8) 


12 


32^ 


10309 


0.241717(6) 


12 


48^ 


16236 


0.241714(3) 


16 


24^ 


15337 


0.176526(6) 


16 


32^ 


18668 


0.176531(3) 


16 


48^ 


19076 


0.1765290(17) 


16 


64^ 


11833 


0.1765302(14) 


20 


24^ 


11484 


0.139295(5) 


20 


32^ 


11634 


0.139283(3) 


20 


48^ 


19814 


0.1392932(13) 


24 


24^ 


15992 


0.115100(3) 


24 


32^ 


20983 


0.1151000(19) 


24 


48^ 


20723 


0.1150986(11) 


24 


64^ 


12101 


0.1151009(9) 


32 


48^ 


20451 


0.0854789(8) 


32 


64^ 


24662 


0.0854815(5) 


32 


96^ 


24875 


0.0854806(3) 


40 


48^ 


20817 


0.0680065(6) 


40 


643 


25442 


0.0680058(4) 


40 


96^ 


25700 


0.06800677(19) 


50 


643 


33448 


0.0541741(3) 


50 


963 


69213 


0.05417428(10) 


50 


1283 


29261 


0.05417418(10) 


50 


3203 


8298 


0.05417406(5) 


64 


963 


25211 


0.04217128(12) 


64 


1283 


35565 


0.04217113(6) 


64 


3203 


7921 


0.04217123(4) 


80 


1283 


34310 


0.03365240(6) 


80 


3203 


8356 


0.03365247(3) 
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SU(4) 


/3 


volume 




(1 - wT^[P])a 


24 


48^ 


79254 


0.2257701(7) 


24 


64^ 


15474 


0.2257703(10) 


32 


48^ 


58752 


0.1651322(6) 


32 


64^ 


16039 


0.1651320(7) 


40 


64^ 


16704 


0.1303851(5) 


40 


96^ 


33574 


0.1303857(2) 


40 


128^ 


32872 


0.13038581(13) 


50 


643 


17257 


0.1033093(4) 


50 


963 


34295 


0.10330876(16) 


50 


1283 


33813 


0.10330879(10) 


58 


963 


34810 


0.08861313(14) 


58 


1283 


33443 


0.08861289(9) 


64 


963 


17342 


0.08007684(17) 


64 


1283 


50908 


0.08007682(6) 


80 


1283 


50664 


0.06372001(5) 


100 


1283 


51510 


0.05076660(4) 



SU(5) 


f3 


volume 




(1 - W^r[P]), 


40 

58 

64 

80 

100 

140 

140 

180 

180 


1283 
1283 
1283 
1283 
1283 
1283 
1603 
1283 
1603 


4667 

8515 

27137 

22616 

20238 

12886 

16049 

12184 

8597 


0.2161236(5) 

0.1447591(2) 

0.13048507(12) 

0.10336875(10) 

0.08208926(8) 

0.05817279(8) 

0.05817267(5) 

0.04505589(6) 

0.04505586(5) 




SU(8) 


/3 


volume 






100 
140 
180 
180 
240 
240 
300 
300 
400 
400 
460 


963 

963 

963 

1283 

963 

1283 

963 

1283 

963 

1283 

1283 


6493 

9789 

7127 

3522 

3755 

3857 

11266 

3831 

18120 

4251 

8656 


0.2285506(5) 

0.1584135(3) 

0.1214678(2) 

0.1214678(2) 

0.0900746(3) 

0.09007497(16) 

0.07160377(11) 

0.07160353(13) 

0.05337892(6) 

0.05337902(8) 

0.04631022(5) 
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